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Abstract 

We consider fractional isoperimetric problems of calculus of variations 
with double integrals via the recent modified Riemann-Liouville approach. 
A necessary optimality condition of Euler-Lagrange type, in the form of a 
multitime fractional PDE, is proved, as well as a sufficient condition and 
fractional natural boundary conditions. 
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1 Introduction 

The calculus of variations was born in 1697 with the solution to the brachis- 
tochrone problem (see, e.g., [IS]). It is a very active research area in the XXI 
century (see, e.g., [TI IT^I^TH^ ). Motivated by the study of several natural 
phenomena in such areas as aerodynamics, economics, medicine, environmental 
engineering, and biology, there has been a recent increase of interest in the study 
of problems of the calculus of variations and optimal control where the cost is a 
multiple integral functional with several independent time variables. The reader 
interested in the area of multitime calculus of variations and multitime optimal 
control is referred to [^[?n[5TH35[l37H55] and references therein. 

Fractional calculus, i.e., the calculus of non-integer order derivatives, has 
its origin also in the 1600s. During three centuries the theory of fractional 
derivatives of real or complex order developed as a pure theoretical field of 
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mathematics, useful only for mathematicians. In the last few decades, however, 
fractional differentiation proved very useful in various fields of applied sciences 
and engineering: physics (classic and quantum mechanics, thermodynamics, 
etc.), chemistry, biology, economics, engineering, signal and image processing, 
and control theory (HHIIIIISSIIMIIM] • 

The calculus of variations and the fractional calculus are connected since 
the XIX century. Indeed, in 1823 Niels Henrik Abel applied the fractional 
calculus in the solution of an integral equation that arises in the formulation 
of the tautochrone problem. This problem, sometimes also called the isochrone 
problem, is that of finding the shape of a frictionless wire lying in a vertical 
plane such that the time of a bead placed on the wire slides to the lowest point 
of the wire in the same time regardless of where the bead is placed. It turns out 
that the cycloid is the isochrone as well as the brachistochrone curve, solving 
simultaneously the brachistochrone problem of the calculus of variations and 
Abel's fractional problem [T]. It is however in the XX century that both areas 
are joined in a unique research field: the fractional calculus of variations. 

The Fractional Calculus of Variations (FCV) was born in 1996-97 with the 
proof, by Riewe, of the Euler-Lagrange fractional differential equations [29ll30] . 
Nowadays, FCV is subject of strong current research - see, e.g., [^Ml fTTll^ . 
The first works on FCV were developed using fractional derivatives in the sense 
of Riemann-Liouville [2]. Later, problems of FCV with Grunwald-Letnikow, 
Caputo, Riesz and Jumarie fractional operators, among others, were consid- 
ered [3l|9l[T2l[20]. The literature on FCV is now vast. However, most results 
refer to the single time case. Results for multitime FCV are scarce, and reduce 
to those in [3lll0ll36] . Here we develop further the theory of multitime fractional 
calculus of variations, by considering fractional isoperimetric problems with two 
independent time variables. Previous results on fractional isoperimetric prob- 
lems are for the single time case only [4j[5] . In our paper we study isoperimetric 
problems for variational functionals with double integrals involving fractional 
partial derivatives. 

The paper is organized as follows. In Section 2 we recall some basic defini- 
tions of multidimensional fractional calculus. Our results are stated and proved 
in Section 3. The main results of the paper include natural boundary condi- 
tions (Theorem 13. Sp and a necessary optimality condition (Theorem 13. 4p that 
becomes sufficient under appropriate convexity assumptions (Theorem 13. 6|) . 

2 Preliminaries 

In this section we fix notations by collecting the definitions of fractional deriva- 
tives and integrals in the modified Riemann-Liouville sense. For more informa- 
tion on the subject we refer the reader to [3l [15HT7l ll9 j . 

Definition 2.1 (The Jumarie fractional derivative [17]). Let / be a continuous 
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function in the interval [a, h] and a E (0, 1). The operator defined by 

X 



r(l — a) dx 

a 

is called the Jumarie fractional derivative of order a. 

Let us consider continuous functions / = /(xi, . . . , x„) defined on 

n 

4=1 

Definition 2.2 (The fractional volume integral 0]). For a £ (0, 1) the fractional 
volume integral of / over the whole domain R is given by 

a\ an 

Definition 2.3 (Fractional partial derivatives [3]). Let xi £ [a^, hi\^i — 1, . . . , n, 
and a G (0, 1). The operator aiD^.[i\ defined by 

Xi 

aiD'^i[i]f{xi, ■■■'^■^'f = Y{1~ a) 'dx^- J ^^^^^^^ /(a^i, • ■ • , x^+i, . . . , x„) 

f ij^l-! • ■ • •! •^i — li '^i+l •,••••) '^n) dt 

is called the zth fractional partial derivative of order a, z = 1, . . . , n. 

Remark 2.1. The Jumarie fractional derivative \15\\17l given by ([T|) can be 

obtained by putting n — 1 in Definition \2.3\ 

X 

aDZWix) = Yl^^Y^ /(^ - " f{a))dt = 

a 

Definition 2.4 (The fractional line integral [3]). Let R = [a, 6] x [c, d]. The 
fractional line integral on dR is defined by 

where 

b 

IaRW = ^ J [fit,c)-fit,d)] [h^tr-'dt 

a 

and 

d 

IorW = « / [/(fe, t) - fia, t)] {d - tr-'dt. 
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3 Main Results 

Let us consider functions u — u{x,y). We assume that the domain of functions 
u contain the rectangle R = [a, b] x [c, d] and are continuous on R. Moreover, 
functions u under our consideration are such that the fractional partial deriva- 
tives aD"[l]u and cDy[2]u are continuous on R, a £ (0, 1). We investigate the 
following fractional problem of the calculus of variations: to minimize a given 
functional 

b d 

J[u{;-)]=a^ J J f{x,y,u,aD^[l]u,,D^[2]u){b~xr-\d-yr-^dydx (2) 

a c 

when subject to an isoperimetric constraint 



b d 




g{x,y,u,aD^[l]u,,D^[2]u){b-xr-\d-yr-^dydx = K (3) 



a c 



and a boundary condition 

u{x,y)\QR = tpix,y). (4) 

We are assuming that ip is some given function, K is a constant, and / and g 
are at least of class of C^. Moreover, we assume that 84 f and d^g have con- 
tinuous fractional partial derivatives a-D"[l]; and 95/ and d^g have continuous 
fractional partial derivatives ci'y [2]. Along the work, we denote by dif and dig 
the standard partial derivatives of / and g with respect to their i-th argument, 
i = 1, . . . ,5. 

Definition 3.1. A continuous function u — u{x,y) that satisfies the given 
isoperimetric constraint ^ and boundary condition is said to be admissible 
for problem ©-g]). 

Remark 3.1. Contrary to the classical setting of the calculus of variations, 
where admissible functions are necessarily differentiate , here we are consider- 
ing our variational problem ([l])-© on the set of continuous curves u (without 
assuming differentiability ofu). Indeed, the modified Riemann-Liouville deriva- 
tives have the advantage of both the standard Riemann-Liouville and Caputo 
fractional derivatives: they are defined for arbitrarily continuous ( not necessar- 
ily differentiable) functions, like the standard Riemann-Liouville ones, and the 
fractional derivative of a constant is equal to zero, as it happens with the Caputo 
derivatives. 

Definition 3.2 (Local minimizer to ©-(HI))- An admissible function u — u{x, y) 
is said to be a local minimizer to problem (O-Q if there exists some 7 > such 
that for all admissible functions ii with ||m — 00 "^7 •^i''^] ~ "^M — 

where 

||u||i,oo := max |w(a;,?/)| + max |ai:'"[l]M(a;,2/)| + mii^ lD^[2]uix,y)\ . 

{x,y)^R {x,y)eR (x,y)£R 
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We make use of the following result proved in [5] : 

Lemma 3.2 (Green's fractional formula [3]). Let h,k, and rj be continuous 
functions whose domains contain R. Then, 

b d 

[h(x, y)aDZ[lMx, y) - k{x, v),D-^[2]7i{x, y)] {b - xT'^d - yT'^dydx 

a c 

b d 

\a — 1 f ^ ^ ,\a — 1 ^ 



[aD^imx, y) D'^[2\k{x, y)] rjix, y){b ~ xT'^d - yT'^dyd. 



X 



+ «! [lUl\{h^)+lU2]{ky^)]. 
Remark 3.3. If rj = on dR in Lemma W^ then 

b d 

[h{x, y)aD^[l]r,ix, y) - k{x, y),D'^[2Ux, y)] (6 - xT'^d - yf-^dydx 



a c 

b d 



[aD^[l]h{x, y) D'^[2]k{x, y)] v{x, y){b - xf'^d - yT-'dydx. 

(5) 

3.1 Necessary Optimality Condition 

The next theorem gives a necessary optimality condition for w to be a solution 
of the fractional isoperimetric problem defined by ([2])- ([4]). 

Theorem 3.4 (Euler-Lagrange fractional optimality condition to ©-(HI)). // 
u is a local minimizer to problem ©-(111), then there exists a nonzero pair of 
constants (Aq, A) such that u satisfies the fractional PDE 

dsH {u} (x, y) -a DZ [l]diH {u} [x, y) D"^ [2]d^H {u) (x, y) ^ (6) 

for all (x, y) £ R, where 

H{x, y, u, V, w, Xq, A) := \of{x, y, u, v, w) + Xg(x, y, u, v, w) 

and, for simplicity of notation, we use the operator {•} defined by 

{u} {x,y) := {x,y,u{x,y),aD2[l]u{x,y),cDy[2]u{x,y),Xo,X) . 

Proof. Let us define the function 

u^{x,y)=u{x,y)+eT](x,y), (7) 

where 77 is such that 77 G C^{R), 

v{x,y)\9R = 0, 
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and e G R. If £ take values sufficiently close to zero, then ([7]) is included into 
the first order neighborhood of u, i.e., there exists 5 > such that Ue 6 Ui{u,S), 
where 

Ui{u, 5) ^^u{x, y) : \\u - < (jj . 

On the other hand. 



uo{x,y) = u. 
Let 



de 



de 



de 



b d 



F{e) = a' / / f{x,y,u,{x,y),aD'^[l]u,{x,y),,D'-y[2]u,{x,y)){b~xr-'{d-yr-'dydx 



and 



b d 



G{e) = a' / / g{x,y,u,[x,y),aD°^[l]u,{x,y),,D^[2]u,{x,y)){b~xr-\d^yr-Hydx 



Define the Lagrange function by 

L{e, Ao, A) = AoF(e) + A (G(e) - K) . 

Then, by the extended Lagrange multiplier rule (see, e.g., [40]), we can choose 
multipliers Aq and A, not both zero, such that 



de ° de 



A 



6 = 



dG 
'dl 



= 0. 



The term ^|^_q is equal to 



e=0 



, y, u„a D^[l]u,,,D'^[2]u,)ib - x^-^d - y^-^' 



(8) 



dydx 



:=0 



b d 



/ / d:if{b-xr~'{d-yr-'dydx 



b d 



+ '^^ J J [ddaD^[l]r, + d5fcD'^[2]7^] {b-xr-\d-yr-^dydx. 

a c 

By ([5]) the last double integral in © may be transformed as follows: 
6 d 

a' J J [d4faD:[l]v + d5fcD'^[2]ii] {b~xr-\d-yr-Hydx 



(9) 



b d 



/ / [aD'i[l]dif+,D'^[2]d^f]Tj{b-xr-^{d-yr~Hydz 
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Hence, 


dF 




de 


e=0 


Similarly, 


dG 




de 


e=0 



b d 



■ j J [dsf -aD^[l]d4f -cD^[2]d5f]r,ib-xr-\d-yr-'dydx. 

a c 

(10) 



b d 



= / / [d^g~aDZ[l\dig-,D'^[2]dr,g\T^{h-xr-\d-yr-Hvdx. 



Substituting and (HU into dH), it results that 

b d 

Ao {d^f -aD%[l]dif -cD'^[2]d^f) 



(11) 



dL{e, Aq, A) 



de 



+ A (935 -a i?^[l]943 -c D'^[2]d^g)\T^{h - - vT-^dydx = 0. 

Finally, since ry = on the fundamental lemma of the calculus of variations 
(see, e.g., [H]) implies that 

d^H {u} (x, y) -a D^[l]diH {u} (x, y) D;;[2]95i? {u} [x, y) = 0. 



□ 



3.2 Natural Boundary Conditions 

In this section we consider problem Q-®, i.e., we consider the case when the 
value of function u — u{x, y) is not preassigned on OR. 

Theorem 3.5 (Fractional natural boundary conditions to ©-([21)). Ifu is a local 
minimizer to problem Q-® ; then u is a solution of the fractional differential 
equation ® . Moreover, it satisfies the following conditions: 

1. diH {u} (a, y) = for all y e [c, d]; 

2. diH {u} (6, y) = for all y e [c, d]; 

3. d^H {u} (x, c) = for all x e [a, 6]; 

4-. d^H {u} (x, d) — for all x G [a, b] . 

Proof. Since in problem ([1])-® no boundary condition is imposed, we do not 
require rj in the proof o Theorem 13.41 to vanish on OR. Therefore, following the 
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proof of Theorem 13. 4[ we obtain 



b d 

2 



{dsH {u} {x, y) +a D'^AmH {u} {x, y) 

+,D'^[2]d^H {u} (x, y)) 7^{b - - yy-^dydx 

+ amji[l][d^H{u}{x,y)ri)+I^MdBH{u}{x,y)in)]^Q, (12) 

where rj is an arbitrary continuous function. In particular, the above equation 
holds for 77 = on dR. If = 0, the second member of the sum in (|12l) 

vanishes and the fundamental lemma of the calculus of variations (see, e.g., [24j) 
implies With this result equation (|12l) takes the form 



diH {u} (6, y)77(6, y){d - y^-'dy - J d^H {u} (a, y)ji{a, y){d - yT'^dy 

c 

b 

dr^H {u} {x, c)f]{x, c){h-x)°'-^dx- I d^H {u} {x, d)r]{x, d){h-x)°'^^ dx = 0. 



(13) 

Let Si — ([a, fe] X c) U ([a, 6] x d) U (6 x [c, d\). Since rj is an arbitrary function, 
we can consider the subclass of functions for which rj{x^y)\g^ — 0. For such 77, 
equation (fTS]) reduces to 

d 



= j diH{u}{a,y)^{a,y){d-yr-^dy. 

C 

By the fundamental lemma of calculus of variations, we obtain that 

diH {u} (a, y) = 

for all y G [c, d] . We prove the other natural boundary conditions in a similar 
way. □ 



3.3 Sufficient Condition 

We now prove a sufBcient condition that ensures existence of global minimum 
under appropriate convexity assumptions. 

Theorem 3.6. Let H{x, y, u, v, w, Aq, A) — Ao/(a;, y, u, v, w) + \g{x, y, u, v, w) 
be a convex function of u, v andw. Ifu{x,y) satisfies (jB]), then for an arbitrary 
admissible function u{-, ■) the following holds: 

J[U{;-)]>J[U{;-)], 

i.e., u{-,-) minimizes ([2l). 
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Proof. Define the following function: 

Kx,y) ■■= u{x,y) - u{x,y). 

Obviously, 

fJ-(.x,y)\gR = 0. 

Since H {u} {x,y) is convex and aD"[l], cDy[2] are linear operators, we obtain 
that 

H{u} {x,y) - H{u} {x,y) 

> {iiix, y) ~ y))d^H {u] (x, y) + {aDZ[l]u{x, y) [l]u{x, y)) O^H {u} (x, y) 

+ {,D'^[2]u{x,y) -,D'^[2]u{x,y)) d^H {u} {x,y) 
= {u{x, y) - u(x, y))d3H {u} (x, y) +a [1] ("(a;, y) - u{x, y)) d^H {u} (x, y) 

+e D^[2] {u{x, y) - u{x, y)) d^H {u} {x, y) 

= tJL{x, y)d3H {u} [x, y) +a y)diH {u} (x, y) D'^ [2]//(x, y)d5H {u} (x, y), 

(14) 

where the Aq and A that appear in |u} (x, y) are constants whose existence is 
assured by Theorem 13.41 Therefore |J 

J[u{.,.)]-J[u(.,.)] 

b d 



a 



f{x, y, u,a DZ [l]u,c [2]u)ib - xT'^d - yT-^ydx 
fix,y, u,a [l]u,, D'^y [2]u) {b - xT-^d - yT-'dydx 



a c 

b d 

+ Xo\a^ J J gix,y,u,a D^[l]u,,D^[2]u){b - xr-\d - yr-'dydx - K 

a c 
b d 

Xo\a' J J g{x,y,u,aD^[l]u,,D^[2]u){b-xr-\d-yr-'dydx-K 

a c 

b d 

■2 /" ^ {H{u}~H{u}){b- x)"-\d-y)''-^dydx. 



^From now on we omit, for brevity, the arguments (x, y). 



9 



Using (dH) and ©, we get 



b d 



a' \ \ {H{u}-H{u}){b-x)°'-'id-y)"-'dydx 



b d 



>aM / fidsH {u}{b-x)"-\d~y)°'~'dyd 



b d 



a" / / {aD'^[l]^Ji^^H{u}+,D^^[2]^Ji^^H{u}){b~xr-\d-yr-^dydz 



a c 
b d 



/ / ^JL^^H{u}{b-x)°'-'{d~y)°'-'dydx 



b d 



■ j j {aDZ [l]diH {u} +, [2]d5H {u}) m(6 -x^-^d- yf-^dydx 



a c 
b d 



J j{d^H{u} D'^[l]diH{u} 

a c 

-cD'^ma^H {u}) ti{b ~ x^-Hd - yr-^dyda 



0. 



Thus, J[U{;-)]>J[U{;-)]. 



□ 



4 Conclusion 

The fractional calculus provides a very useful framework to deal with nonlocal 
dynamics: if one wants to include memory effects, i.e., the influence of the past 
on the behaviour of the system at present time, then one may use fractional 
derivatives. The proof of fractional Euler-Lagrange equations is a subject of 
strong current study because of its numerous applications. However, while the 
single time case is well developed, the multitime fractional variational theory 
is in its childhood, and much remains to be done. In this work we consider a 
new class of multitime fractional functionals of the calculus of variations subject 
to isoperimetric constraints. We prove both necessary and sufficient optimality 
conditions via the modified Riemann-Liouville approach. 
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